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Abstract. It is shown that the new formula for the held theory Poisson brackets arise naturally 
in the extension of the formal variational calculus incorporating divergences. The linear spaces 
of local functionals, evolutionary vector fields, functional forms, multi-vectors and differential 
operators become graded with respect to divergences. The bilinear operations, such as the action 
of vector fields on functionals, the commutator of vector fields, the interior product of forms and 
vectors and the Schouten-Nijenhuis bracket are compatible with the grading. A definition of the 
adjoint graded operator is proposed and antisymmetric operators are constructed with the help 
of boundary terms. The fulfilment of the Jacobi identity for the new Poisson brackets is shown 
to be equivalent to vanishing of the Schouten-Nijenhuis bracket of the Poisson bivector with 
itself. It is demonstrated, as an example, that the second structure of the Korteweg-de Vries 
equation is not Hamiltonian with respect to the new brackets until special boundary conditions 
are prescribed. 

1. Introduction 

The Hamiltonian formalism of classical mechanics may serve as the ideal model illustrating 
the harmony of physics and mathematics. During the last 20 years it was realized that a 
number of its mathematical constructions, for example, the Schouten-Nijenhuis bracket [1], 
could be extrapolated to field theory[2, 3]. This made the search for new nonlinear integrable 
models much easier. Even more general constructions uniting the Schouten-Nijenhuis bracket 
with the Frolicher-Nijenhuis bracket were considered by A. Vinogradov [4]. 

But the methods mentioned above (and usually called the formal variational calculus[5]) 
have some restrictions arising from boundary conditions which should allow free integration 
by parts. As a rule, the periodic boundary conditions or the rapid decay of fields at spa- 
tial infinity are necessary. Of course, these are not all the physically interesting cases. For 
example, the Coulomb potential in electrodynamics does not tend to zero rapidly enough. 
The similar behaviour is typical for Yang-Mills and gravitation fields. Nontrivial boundary 
problems arise also in the motion of material continua. 
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We consider the Hamiltonian treatment of nontrivial boundary problems as an interesting 
field of research where there is some place for new approaches and results. The field theory 
Poisson brackets which fulfil the Jacobi identity under arbitrary boundary conditions have 
been proposed in [6]. Here we extend the formal variational calculus to the most general case 
when no one boundary term arising in integration by parts can be discarded. We hope to 
present some physical applications of the methods developed here in the future. Let us also 
say that the interest to the role of divergences in field theory is vivid now as can be seen 
from related papers [7, 8, 9]. 

As an example, illustrating the nonstandard nature of the problems to be considered, let 
us remind the history of the long discussion on the role of surface integrals in the canonical 
formulation of General Relativity. During about 15 years Arnowitt, Deser and Misner[10], 
Diracfll], Higgs[12], Schwinger[13], DeWitt[14], Regge and Teitelboim[15] were involved in it. 
The solution obtained in the last work[15] serves as the paradigm for the treatment of similar 
problems up to now. It has been proposed to work with the special class of the so-called 
differentiable or admissible functionals. These functionals are defined by the requirement 
that their variation should not have surface contributions under the prescribed boundary 
conditions. The Poisson brackets for these functionals are the standard ones, i.e., they are 
just the same as given in the formal variational calculus 



but now nonzero surface contributions are allowed. 

Here a natural question to ask is: do these brackets fulfil the standard axiomatic require- 
ments, i.e., the Jacobi identity and the existence of the Poisson algebra on this space of 
admissible functionals? For the infinite domain of integration and the asymptotic boundary 
conditions the affirmative result for the second requirement was obtained by Brown and 
Henneaux[16]. The first requirement was partially analysed by us and in the case treated 
above the answer is also positive. 

It is more difficult to study the finite domain case. Let us take as a second example the 
motion of a fluid or plasma. It was shown by Lewis, Marsden, Montgomery and Ratiu[17] 
that the Jacobi identity for the standard Poisson bracket can be violated even in the case 
of fixed boundary, and so the Poisson brackets should be modified by surface terms. In the 
free boundary case it turns out to be natural to extend the space of admissible functionals 
so that their variation could include nonzero surface contributions. But according to [17] the 
presence of nonzero term with 5q A in the boundary integral requires the absence of the 
corresponding term with Spa and vice versa. A new formula for Poisson brackets arises as 
a result of a generalization of the variational derivative which is now allowed to contain a 
surface contribution 




5q A (x) 5pa(x) 5q A (x) Spa(x) 



SF 5G 5G SF 



d n x, 



n 



x + 




m 



X + 



J ~5pA~ SpA ^ dndS ' 
an 



n 



Unfortunately, it is not quite clear whether the Poisson bracket of the two functionals, 
admissible in this new sense, will be admissible functional itself. 
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As a third example, we would like to attract reader's attention to consequences of the non- 
commutativity of the standard variational derivatives, i.e., the Euler-Lagrange derivatives. 
This point was discussed formally in publications by I. Anderson[19, 20] and Aldersley[21]. 
We faced with the problem independently, in the course of studying surface terms in the 
Poisson algebra of the Ashtekar formalism of canonical gravity [18]. There it was found that 
transformations of the type 

q A (x) -> q A (x), p A (x) -> /Ml.'') + 5/1 



Sq A (x) 



were noncanonical ones if surface terms were not ignored. Tracing the connection of this 
calculation with the standard calculations with <5-functions[22], we have found that the cor- 
respondence could be restored by introduction of 6q — the characteristic function of the 
domain Q 



a ( \ f 1 if x e 
0n(aO = { 



Then the standard relations 



otherwise. 



+ — )5(x,y)=0, 



dx l dy l 



d d \ d6 

On(x)— + 9 n (y)— 5(x,y) = -—^8{x,y), 



should be modified as 



dx % dy l J dx 

where the usually discarded surface term is preserved. 

All the above examples lead us to the necessity to extend the formal variational calculus 
onto total divergences. This extension consists in the introduction of a new grading for the 
linear spaces of local functionals, vector fields, functional forms, multi- vectors and differential 
operators. To come back to the standard case one should put 6q{x) = 1 in W 1 . 

The extension of the formal variational calculus naturally incorporates the new definition 
of local functionals (not modulo divergences) and of their differentials (as a full variation, 
not fixed on the boundary). The Poisson bivectors are now defined in a more general way 
and they can include boundary contributions. The definition of pairing (interior product) 
is also revised and now the trace of two differential operators is used for it, so the pairing 
becomes compatible with the grading. The Poisson bracket found more or less heuristically 
in [6] arises now on the base of the geometrical constructions as 

{F,G} = dGjdFj^, 

where \& is the Poisson bivector. 

We show here that the Jacobi identity for the new Poisson brackets can be verified 
without the long calculations of binomial sums used in [6]. Its fulfilment is equivalent to 
the vanishing of the Schouten-Nijenhuis bracket of the Poisson bivector with itself. And in 
its turn this condition can be easily tested along with the procedure proposed by 01ver[2] 
with a minimal modification of it. More attention than in [6] is paid here to nonultralocal 
Hamiltonian operators with nonconstant coefficients, because now the technical obstacles are 
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removed. It turns out to be that not all operators which are Hamiltonian with respect to 
the standard brackets remain Hamiltonian in relation to the new brackets. For example, the 
second structure of the Korteweg-de Vries equation is not automatically Hamiltonian with 
respect to the new formalism. In this respect it strongly differs from the first KdV structure. 

The content of this work is as follows. In Section 2 we introduce the grading for the local 
functionals, and the evolutionary vector fields. In Section 3 the differential, the functional 
m-forms, the interior product of vector fields onto m-forms and Lie derivative are defined. 
Section 4 deals with graded differential operators and their adjoints. In Section 5 we discuss 
multi-vectors and the Schouten-Nijenhuis bracket. It is shown how 1-vectors and evolutionary 
vector field spaces are isomorphic. Section 6 contains the general definition of the Poisson 
bracket, its realization in this formalism, the definition of Hamiltonian vector fields and 
the theorem on connection between the commutator of two Hamiltonian vector fields and 
the Poisson bracket of corresponding Hamiltonians. All constructions are illustrated by an 
example: the first Hamiltonian structure of the Korteweg-de Vries equation. The proof of 
Jacobi identity is postponed until Section 7. This new proof is applicable for all local Poisson 
brackets and so supersedes the proofs given earlier in [6]. In the same time it is easy to 
compare this proofs because notations are the same. At last, in Section 8 we consider two 
examples of the non-ultralocal Poisson brackets with field dependent coefficients (this class of 
brackets we were unable to study by the methods of [6]). We show that the second structure 
of Korteweg-de Vries equation is not Hamiltonian if surface terms are not ignored, whereas 
the 2-dimensional flow of the ideal fluid is described by Hamiltonian structure. This points 
to nontrivial character of the generalization of the formal variational calculus suggested here. 
A short summary is given in Conclusion and the points remaining unclear are mentioned. 
We hope to continue this research by the detailed treatment of the boundary conditions role 
and applications to some physical problems. 

As a rule, we use below the same notations as in [6] except a change of the notation for 
the Frechet derivative from Df to /' and omitting the sign of summation according to the 
Einstein rule. We find convenient to represent integrals over finite domain Q through integrals 
over the infinite space R n by inserting into all integrands the characteristic function 9q. Then 
the formalism seems closer to the standard formal variational calculus where local functionals 
and functional forms are defined modulo divergences. But the formal divergences that we 
discard here are integrated to zero under arbitrary conditions on the boundary of the finite 
domain, whereas real divergences are incorporated into graded structures. All the operations 
introduced below are compatible both with discarding formal total divergences (if one object 
is a formal divergence than the result of operation is also formal divergence) and with the 
grading (i.e., the same is valid for real divergences). Extension of the space of differential 
operators by admitting their grading permits to use the concept of adjoint operator. So, 
antisymmetric operators can now be constructed and the Poisson bracket formulas become 
more compact, than in [6] though their content is the same. Nevertheless, in the proof of the 
Jacobi identity we prefer to use the old notations to make easier the comparison with the 
not so general proofs of [6]. 

Henceforth we consider the space M. n and the multi-index notations J = (ji, ...,j n ) where 
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ji > 

AJ) _ 



A &> 1 x 1 ...&> n x n '' 
The Frechet derivative is defined as 



1^1 =jl + ••• + jr. 



df 



where 

n - — -+- 6^ J+i) _ n r - n jl n jn n° - i - n 

Binomial coefficients for multi-indices are 

K) = Gi)"'Gn)' 

A = |j!/W-fc)0 if < A: < j; 
fc/ lo otherwise. 

With the help of them we introduce the so-called higher Eulerian operators [2, 21, 23] 

EiU) = (-D^QD K ^. (2) 



2. Local functionals and evolutionary vector fields 



Let us start with notions from the theory of graded spaces as they are given in Ref. [3]. A 
grading in linear space L is a decomposition of it into direct sum of subspaces, with a special 
value of some function p (grading function) assigned to all the elements of any subspace. 

Below the function p takes its values in the set of all positive multi- indices J = (ji, . . . , j n ) 
and so, 

oo 

£ = 0£ <J> . 

J=0 

Elements of each subspace are called homogeneous. 

A bilinear operation x, y h- > x o y, defined on L, is said to be compatible with the grading 
if the product of any homogeneous elements is also homogeneous, and if 

p(xoy) =p(x) +p(y). 

Now let us turn to the concrete structures. 

There are two ways to write a local functional: as the integral of a smooth function 
/(°) [(^a\ x )) °f fields and their derivatives up to some finite order over the prescribed 
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domain f2 in W 1 , or as the integral over all the space W 1 but with the characteristic function 
of the domain Oq, included into the integrand 

F = J f {0) (<pf\x)) d n x = J 9 n f^d n x. (3) 

As in [6] , let us denote the space of local functionals as A. Here we shall call the expression 
given above the canonical form of a local functional. We formally extend that definition by 
allowing local functionals to be written as follows 



Dj6 n f {J) {<f>f\x))d n x = j 6^f^d n x = j fd n x, 



where only a finite number of terms is allowed. Here and below we simplify the notation for 
derivatives of 9 and remove fl All the integrals without the domain of integration shown 
explicitly are integrals over W 1 , below we shall omit d n x. Of course, any functional can be 
transformed to the above form (3), exclusively used in [6], through integration by parts 



where 

;<o> = { -if\ Djf w. 

Evidently, the formal integration by parts over infinite space W 1 changes the grading. It 
will be clear below that the general situation is the following: from one side we have the 
compatibility of all the bilinear operations with the grading and from the other side - 
with the formal integration by parts. So, basic objects (local functionals etc.) are defined 
as equivalence classes modulo formal divergences (i.e., divergences of expressions containing 
^-factors) and the unique decomposition into the homogeneous subspaces with the fixed 
grading function can be made only for representatives of these classes. 
We call expressions of the form 



the evolutionary vector fields. The action of the evolutionary vector field on a local functional 
is given by the expression 

*F = J 6^D K ^^ = J 9^ J V J) f {I) - / i>f- (4) 

It is a straightforward calculation to check that this operation is compatible with the formal 
integration by parts, i.e. 

4Df = Dtyf), 
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as it is in the standard formal variational calculus. This relation is, of course, valid for 
integrands. 

It is easy to see that the evolutionary vector field with coefficients 



^ - Dl ^ b - DlXb 



'3 ^r'B 

can be considered as the commutator of the evolutionary vector fields S and A 



M>F=[S,A]F = J(t(\f)-\(tfj), 



with the Jacobi identity fulfilled for the commutator operation, and so these vector fields 
form a Lie algebra. 

Let us comment upon the representation of the evolutionary vector fields as integrals, 
which is different from the traditional notations. 

The formal variational calculus [5] operates with the local functionals which are expressed 
by single integrals of functions of the specified class, for example, infinitely differentiable 
ones. The functional forms and multi-vectors are expressed by similar integrals. The pairing 
of two such objects gives us a single integral again. 

At the same time, another notations are widespread, especially in physical literature, 
which use 5-function and its derivatives. Then a result of the pairing of two single integrals 
is understood as a double integral. But as this double integral contains the 5-function, it 
always can be converted into the single one. 

This convertation of a double integral into the single one with the help of (^-function is 
trivial when no boundary terms could arise. The subject of this work is just a study of the 
opposite case. The new rule is necessary here and it have been proposed in [6] as Rule 4.2 

f(x)g(y)Df ) D^S(x,y) = J D K fD j9 . (5) 




n n 



In this article we give really a new and equivalent form of these rules which help to avoid 
the usage of double integrals at all. 

The concept of vector field appeared initially in the course of studying the evolutionary 
differential equations and their symmetries. In the formal variational calculus [5] functionals 
are, in fact, replaced by equivalence classes of functions, and so the action of evolutionary 
vector fields onto local functionals is replaced by their action on functions 

df 



ipf = D K ip A - 



However, to represent functionals by integrals and to require that the result of the action of 
an evolutionary vector field onto a local functional is a local functional, i.e. an integral, it is 
absolutely natural to represent the evolutionary vector fields also as integrals 



d_ 
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in combination with the standard rule 



d<pf\x) 



S(x,y)5 AB 8 JK . 



(6) 



Another argument supporting our notations is the equivalence between evolutionary vec- 
tor fields and 1-vectors, which is demonstrated for the standard formal variational calculus 
in book [2] and also for the graded case in Section 5 of this article. 1-vectors as a partial case 
of multi- vectors are always written as integrals. 

Apart from the notational revision we would like to mention a new feature in our treat- 
ment of the vector fields: now they are not derivations when applied to standard functions, 
but only to the graded ones. Of course, in the traditional approach the vector fields are not 
derivations when applied to functionals, because their multiplication is not defined. But these 
vector fields, traditionally written without the integral sign, are derivations of functions. This 
property is partially lost here. It can be restored formally if we consider integrands containing 
9 as functions and take a relation 



as a definition of their multiplication. 

In this context, formula (4), introduced as a definition, can be interpreted also as a 
consequence of the standard relation (6) and a new definition (7). 

Therefore, it is evident that our "rule for multiplication of distributions" taken from [6], 
i.e., equation (7) is nothing more than another way to define the pairing compatible with 
the introduced grading. 

At last, let us mention the possibility to use other notations in this formalism. It is, of 
course, possible to avoid ^-functions and to use integrals over the domain Q only. Then any 
local functional can be given as 



with analogous rewriting of the other objects. Correspondingly, equation (4) will be written 



Dj6 x Dj6 = D I+J 6. 



(7) 




where 



f ((J)) = (_l)l^>, 



as 




3. Differentials and functional forms 



The differential of a local functional is simply the first variation of it 
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here and below 8<j)^ = D K 5(f) A . It can also be expressed through the Frechet derivative (1) 
or through the higher Eulerian operators (2) 

dF = J eWfW'w) = J e^D K (E«{fV))5<t> A ) . 

This differential is a special example of functional 1-form. A general functional 1-form can 
be written as 



S = j flW^W = J e (J) a (J) = J a . 



Of course, the coefficients o^ AK are not unique since we can make formal integration by parts. 
Let us call the following expression the canonical form of a functional 1-form 



a J) Ha- 



Analogously, we can define functional m-forms as integrals or equivalence classes modulo 
formal divergences of vertical m-forms 

= = s / A.-.A = / *V» = / o. 

Define the pairing (or the interior product) of an evolutionary vector field and 1-form as 

m = s J s = J e^a A \D K ^ = f ev+W(t<n) = J a®. (8) 

The interior product of an evolutionary vector field and a functional m-form will be given 
as follows 

s j s = JL(-i)« / e-'+Mi- a ,k„Dk,(^' a... 

...A^' ) A^'>A...A^~». 
Then a value of the m-form on the m evolutionary vector fields will be defined by the formula 

. . . , n m ) = ^ m _| . . . E. 

It can be checked by straightforward calculation that 

( j Dct)(6,...,U = ^(^(6,...,U)- 

The differential of the m-form which is given as 

dE = _l j 6 {j) d °MK^K m 5 $<) A ^) A ... A ^(£0 = y e (j) d(T <j> = y d(7> 

satisfies standard properties 

d 2 = 
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and 

dS(^!, . . . , c m+1 ) = ^(— l) l+1 ^jS(r, 1 , . . . , Cj, . . . , c m+1 )+ 

i 

+ l) t+J S([^, ci, . . . , Cj, . . . , Cj, . . . , c m+ i). 

The Lie derivative of a functional form X along the evolutionary vector field S can be 
introduced by the standard formula 



L H S = Hjd£ + d(Hj£). 



4. Graded differential operators and their ad joints 
We call linear matrix differential operators of the form 

Iab = D N 

the graded differential operators. 

Let us call the linear differential operator I* the adjoint to / if for an arbitrary set of 
smooth functions /a, 9 a 

J IaIabOb = J QaIabIb- 
For coefficients of the adjoint operator we can derive the expression 

r}S M = ( _ 1)|Jt| Q (K- Ly K _ L _ M ,u-m (9) 

It is easy to check that the relation 

I A B(x)5(x,y) = I* BA (y)5(x,y) 
follows from Rule 4.2 of [6]. For example, we have 

In one of our previous publications [18] we tried to connect the appearance of surface terms 
in Poisson brackets and the standard manipulations with the 5-function. The ansatz used 
there for the above simplest example coincided with (10) up to the sign. The reason for this 
difference laid in the other choice made there instead of Rule 4.2 of [6]. That ansatz lead us to 
the standard Poisson brackets which are not appropriate for nontrivial boundary problems. 
Operators satisfying the relation 

r = -i 
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will be called the antisymmetric ones. With the help of them it is possible to express 2-forms 
(and also 2-vectors to be defined below) in the canonical form 




It is clear that we can consider representations of functional forms as decompositions over 
the basis derived as a tensor product of 5<pA, with the totally antisymmetric multilinear 
operators 

as coefficients of these decompositions. 



5. MULTI- VECTORS, MIXED TENSORS AND SCHOUTEN-NlJENHUIS BRACKET 

Let us introduce dual basis to \54>a) by the relation 

,5<j) A (x) ) = S AB S(x,y) (11) 



and construct by means of the tensor product a basis 

5 5 



Hb 1 (v) S<f>B 2 (y) HbSv)' 

Then by using totally antisymmetric multilinear operators described in the previous Section 
we can define the functional m-vectors (or multi-vectors) 

* = ^i/^ _^A...A^ = /^> 

Here a natural question on the relation between evolutionary vector fields and 1-vectors 
arises. Evidently, evolutionary vector fields lose their form when being integrated by parts 
whereas 1-vectors preserve it. Let us make a partial integration in the expression of a general 
evolutionary vector field 



by removing D K from C, A J \ then we get 



s = |ei J> ^ (J+L) (-i) |x| ( 



A 



K\ ^ d 



(K)- 



It is easy to see that by using (7), i.e., Rule 5.4 from [6], in the backward direction we can 
write 



Ha 
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where the higher Eulerian operators (2) and the full variational derivative (Definition 5.1 of 

[6]) 

are consequently used. Therefore, we arrive at the following Statement. 

Statement 5.1 There is a one-to-one correspondence between the evolutionary vector 
fields and the functional 1-vectors. The coefficients of 1-vector in the canonical form ^ are 
equal to the characteristics of the evolutionary vector field. 

It is not difficult to show that we can deduce the pairing (interior product) of 1-forms 
and 1-vectors and this pairing preserves the identification. Really, the definition of the dual 
basis (11) and (7), i.e., Rule 5.4 of [6], permits us to derive that 

S(2) = SJS = J J eU(x)eW(y)*V(x)Zjp L (y) (d l1 ^,D k 6M*)) = 

eV+»D L a<&D K t$ i = J eV+Wtf*) = J (7(0 = J ^+ J )Tr(^>^>), 

and when 1-vector is in the canonical form (only L = term is nonzero) this result coincides 
with Eq.(8). 

This formula for the pairing will be exploited below also for interior product of 1-vectors 
and m-forms or 1-forms and m- vectors. Its importance comes from the fact that it is invariant 
under the formal partial integration both in forms and in vectors, i.e., 

(D(7)(0=D((7(0)=(7(D(0). 

Evidently, it is the trace construction for convolution of differential operators (as coefficients 
of tensor objects in the proposed basis) that guarantees this invariance. 

The interior product of 1-vector onto m-form and, analogously, of 1-form onto m-vector 
is defined as 

hjs = ^-ir/^o^fa w*X''a... 

... AW£r>A^>A...AW£">) -(-1)< H,) / Jl». (12) 

Then we also can define the value of m-form on m 1-vectors (or, analogously, m-vector on 
m 1-forms) 

■ ■ • , = S m J . . . Sx J E = J 0(^+-+^)Tr (ct^ ' • ' d /m> ) , 

where each entry of multilinear operator a acts only on the corresponding £, whereas each 
derivation of the operator £ acts on the product of a and all the rest of £'s. 
It is possible to define the differential of m-vector 

M = ~, [ 9^ ^ K ^- S^D Kl ^- A ... A D Km J-, 
ml J S< t>M 5<t>A m 
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as an example of a mixed (™) object. Evidently, d 2 ^ = 0. 

With the help of the previous constructions we can define the Schouten-Nijenhuis bracket 

[S,*] SJV = dSj* + (-l) M d*jS 

for two multi-vectors of orders p and q. The result of this operation is p + q — 1-vector 
and it is analogous to the Schouten-Nijenhuis bracket in tensor analysis [1]. Its use in the 
formal variational calculus is described in Refs.[2, 3]. However, in cited references this bracket 
is usually defined for operators. We can recommend Ref.[24] as an interesting source for 
the treatment of the Schouten-Nijenhuis bracket of multi-vectors. Our construction of this 
bracket guarantees a compatibility with the equivalence modulo divergences 

m^} SN = m^} S N = ^D^ SN . 

Statement 5.2 The Schouten-Nijenhuis bracket of functional 1 -vectors up to a sign coincides 
with the commutator of the corresponding evolutionary vector fields. 

Proof Let us take the two 1-vectors in canonical form without loss of generality 




and compute 

[~,*] 5iV = d~J*-d*J~. 

We have 



and 




Therefore, we obtain 

and the proof is completed. 

Statement 5.3 (Olver's Lemma [2]) The Schouten-Nijenhuis bracket for two bivectors 
can be expressed in the form 

[a, *\ SN = -\j £ a i\m H-\f e a k\k) a e, as) 

where the two differential operators I, K are the coefficients of the bivectors in their canonical 
form. 
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Proof. Let us consider the Schouten-Nijenhuis bracket for the two bivectors and without 
loss of generality take them in the canonical form 

A = lJe iL) UAi^ N D N £ B , 

where £a = 5 /5(f) a and operators I , K are antisymmetric. Then we have 

* J d(f) y E 

and 

-\ J e^D P {j^U A D N ^ A Dj (ZcKg2 P ) . 
Now let us make integration by parts in the second term 

dA J * = -1 f A (l l A T)' (K' M> () A D N £, B - 

_1 J <L + M +QH _ lfl Q a^_ u A Dn(b a Dj+pQ / fcJf (M>P) 

At last we change the order of multipliers under wedge product in the second term, make a 
replacement M — > M — Q and organize the whole expression in the form 

dA J * = -\ J e^u a (j<g")' o (4?&+ 

Having in mind the definition of adjoint operator (9) we can represent the final result of the 
calculation as follows, 

[A,*w = -\J o (L+M kA ((i^)'(K^o + (K {M W {L) 0)^Z, 

thus supporting in this extended formulation the method, proposed in Ref. [2] for testing the 
Jacobi identity (see Section 7). For the general procedure of testing Hamiltonian properties 
see also [26]. 
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6. POISSON BRACKETS AND HAMILTONIAN VECTOR FIELDS 



Let us call the bivector 

. 1 f 5 * 5 




formed with the help of the graded antisymmetric differential operator 



Iab = 0^1$? D. 



the Poisson bivector if 




{F, G} = *(dF, dC7) = dG J dF J * 



the Poisson bracket of these functional. 

The explicit form of the Poisson bracket can easily be obtained. It depends on the explicit 
form of the functional differential, which can be changed by the formal partial integration. Of 
course, all the possible forms are equivalent. Taking the extreme cases we get an expression 
through Frechet derivatives 



Theorem 6.1 The Poisson bracket defined above satisfies the standard requirements of 
the bilinearity, antisymmetry and closeness on the space of local Junctionals A, i.e., Definition 



Proof. 1) From the previous formulas (14), (15) it is clear that {F, G} is a local functional, 
2) antisymmetry of {F, G} is evident and 3) equivalence of the Jacobi identity to the Poisson 
bivector property will be proved in Section 7. 

The result of interior product of the differential of a local functional H on the Poisson 
bivector (up to the sign) will be called the Hamiltonian vector field (or the Hamiltonian 




(14) 



or through higher Eulerian operators (2) 




(15) 



2.3 of [6]. 



1-vector) 



IdH = -dHjy 



corresponding to the Hamiltonian H. Evidently, the standard relations take place 



{F, H} = dF(IdH) = (IdH)F. 
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Theorem 6.2 The Hamiltonian vector field corresponding to the Poisson bracket of the 
junctionals F and H coincides up to the sign with commutator of the Hamiltonian vector 
fields corresponding to these functionals. 

Proof. Consider a value of the commutator of Hamiltonian vector fields IdF and IdH on 
the arbitrary functional G 

[IdF,idH]G = IdF(IdH(G)) - IdH(IdF(G)) = 

= IdF({G, H}) - IdH({G, F}) = {{G, H}, F} - {{G, F}, H} = 

= -{G,{F,H}} = -Id{F,H}(G), 

where we have used the Jacobi identity and antisymmetry of Poisson bracket. Due to the 
arbitrariness of G the proof is completed. 

Example 6.3 

Let us consider the first structure 

{u(x),u(y)} = \{D X - D y )5(x,y) 

of the Korteweg-de Vries equation (Example 7.6 of Ref. [2]) 

Uf u xxx -\- uu x . 

Construct the adjoint graded operator to 6D according to Eq.(9) 

(9D)* = -9D - DO 

and the antisymmetric operator is 

I = ±.(0D- (6D)*^ =9D + ^D9. 
The Poisson bivector has a form 

* = i fed* Dpi. 

2 J \Su 5uJ 

The differential of a local functional H (for simplicity it is written in canonical 




form) is equal to 

dH = Jeti{5u) = J 9 {k \-l) k E k (h)5u, 

where the Frechet derivative or higher Eulerian operators can be used. Therefore, the Hamil- 
tonian vector field generated by H is 

^ = ^J* = ~/»(V(^)-«V(A)V 
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or 

5 



_I J e^ k \-l) k (E\h)D - DE\h)) 



5u 
or also 

-\ J 6> (fc) (-l)*A {E\h)D - DE\h)) JL. 
The value of this vector field on another functional F coincides with the Poisson bracket 
-dFjdi/J^ = {F,H} = ^ J 0( k + l \-l) k + l (E k (f)DE l (h) - E k {h)DE\f)) . 



7. Proof of Jacobi identity 

In this section we will prove that the Jacobi identity for the Poisson bracket is fulfilled if 
and only if the Schouten-Nijenhuis bracket of the corresponding Poisson bivector with itself 
is equal to zero. This should complete the proof of Theorem 6.1. 

Let us use one of the possible forms of the Poisson brackets given in Appendix of [6] 

{F,G} = \J ^Tr (/'(/ <J V) " g\i {J) f)) , 

where the differential operator I is not supposed to be antisymmetric for the easier compar- 
ison of this proof with that given in [6] . We remind that in less condensed notations 



(/W)) 



Tl ' ( '" l '"' ,) ^ [ m){ K l) D ^ D ^- »^VW, Jjy 



(in Appendix of [6] the indices M and L in the binomial coefficients of the same formula are 
unfortunately given in the opposite order). 
We will estimate the bracket 

{{F,G},H}= l -J 9^TT[{f,g}\i^h')-h'(i^{f,gy)), 
where {/, g} denotes the integrand of {F, G}. Since Frechet derivative is a derivation we have 

{/, g}' = ^<*>Tt(r(/<*V, •) + fi' {K) {-)g' + g"{fi {K \ •)-(/<- g)) 



and 



Tr ({/, g}'ih') = \ (nig', Ih>) + f'P(ih')g' + g"(f'I, Ih>) -(/<-> g)) . 



Let us explain that /" denotes the second Frechet derivative, i.e., the symmetric bilinear 
operator arising in calculation of the second variation of the local functional F (in the 
canonical form): 
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When we put into entries of /" operators under the trace sign it should be understood that 
these operators act on everything except their own coefficients, for example, 

Trir W j h ^ = ^( L Z p ]( M V M - s u 



(nh'Jh') 



PJ\ Q J\S J\ T 



d 2 f ^ [ n j dg \ ^ ( n ; dh 

i — 

and the expression remains symmetric under permutation of its entries 



xZW-p-q-s-t {J) 7Tak)Dj + t DpIac—m D k+q A^d— ^ 



tv (r (V, /V)) = tc (r (/v, /V)) 



When the operator J stands to the right from the operator of Frechet derivative /' as in 
expression 

TV (g"(Ih', /'/)), 

it acts on everything except /'. At last, for Frechet derivative of the operator we have 



Making similar calculations we get 

Tr (h'I{f, g}') = X -Tr (f"(h'I, Ig>) + fi'(til)g' + g"(f'I, til) - (/ - «?)) 
and therefore 

{{F,G},H} = -^J e^ J+K ^T^f'(i^g'J^ti) - f"(tii( J \iWg')- 

-f"(i^h', g 'i^) + f"( g 'iv\h'iW) + fi'^(i^ti - tii^) g > - (/ - A 



Just the first four terms, apart from the fifth containing Frechet derivative of the operator 
I, were present in our proof for nonultralocal case given in [6] (only terms with zero grading 
were allowed for / there). After cyclic permutation of F, G, H all terms with the symmetric 
operator of the second Frechet derivative are mutually cancelled and 

{{F,G},H} + c.p. = ^ e^ J+K ^T^fP^(i^ti - h'iW)g'- 

-g'i'( J \i^ti-tii^)f' + c. P )j, 

where cyclic permutations of F, G, and H are abbreviated to c.p.. When operator / is given 
in explicitly antisymmetric form all the four terms are equal. Taking into account Olver's 
Lemma (13) we get 

{{F, G}, H} + c.p. = -[/, I] SN (dF, dG, dH), 

so finishing the proof. 
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8. Examples of nonultralocal operators 

The second structure of the Korteweg-de Vries equation may serve as a counter-example to 
the hypothesis [25] that all operators which are Hamiltonian with respect to the standard 
Poisson brackets should also be Hamiltonian in the new brackets. 

Example 8.1 

Let us start with the standard expression (Example 7.6 of Ref. [2]) 
...... / d 3 2d 1 du\ . 

{^)My)}=[-^ + - 3 u Tx + - 3Tx ) S(x,y) 

and construct the adjoint operator to 

K = 9 (d 3 + ^uD + ^Dv)j , 

which is 



K* = -9 {D 3 + -uD + -Duj - D9 y$D 2 + -uj - 3D 2 9D - D 3 9. 
Then the antisymmetric operator 

I = \{K- K*) = 9^D 3 + ^uD + ^Du^j + D9 (^D 2 + ^ + ^D 2 9D + l -D 3 9 
can be used for forming the bivector 

* = i/«A/{, 

where 5/5u = £. This bivector has a form 

* = \ J (ot A D£ + 1d6£ A D 2 i + {^D 2 9 + 2 -9u)t A D^j . 

Then evaluating the Schouten-Nijenhuis bracket for the bivector with the help of Statement 
5.3 

and integrating the first term by parts we get 

Therefore, instead of the Jacobi identity we have 

{{F, G}, H} + c.p. = - l -J D i+j+k+1 (E i (f)DE j (g)D 2 E k (h) + c.p.) dx. 



n 
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So, the second structure of KdV equation can be Hamiltonian only under special boundary 
conditions. 

Example 8.2 

Now let consider another example which is also nonultralocal, but the operator remains 
to be Hamiltonian in the new brackets independently of boundary conditions. The Euler 
equations for the flow of ideal fluid can be written [2] in Hamiltonian form as follows (Example 
7.10 of Ref. [2]) 

du _ 5H 
dt 5uo ' 

where 

H— I -\\i\ 2 d 2 x, oj = V x u. 



T 

Let us limit our consideration by the 2-dimensional case when uj has only one component u 
and 

V = U! X Dy — U!yD x , 

where Ui = DiU, i = (x,y). We can construct the antisymmetric operator 

i= l -(6V- (6V)*} = 0{u x D y - u y D x ) + l -{D y 6u x - D x 9u y ), 
and then the bivector 

where £ = 5/5uj. Statement 5.3 gives us 

[*. ^]sn = J (o MZ A Cxy A i y - i A iyy A + u y (Z A i X y A - £ A £ xx A Q) + 

+ (D y 6u x - D x 6u y ) £ A A ^ 
and after integration by parts the expression can be reduced to zero. 



9. Conclusion 



We have shown that there is an extension of the standard formal variational calculus which 
incorporates divergences without any specification of boundary conditions. It should be im- 
portant to understand relations of this formalism to the constructions of the variational 
bicomplex [27]. It seems also rather interesting to study if some physically relevant algebras 
can be realized with the help of the new Poisson brackets as algebras of local functionals. 
One such example is considered in [28]. 
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